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Longhi et al. [1] have used the D(Ni) vs. D(Mg) 
parameterizations of Jones [2, 3] in attempting to ex- 
plain the Ni systematic s of lunar differentiation. A key 
element of the Jones parameterization and the Longhi 
et al. models is that, at very high temperatures, Ni may 
become incompatible in olivine. Unfortunately, there 
is no actual experimental evidence that this is ever the 
case [1]. To date, all experiments designed to demon- 
strate such incompatibility have failed. 

Here I will investigate the thermodynamic founda- 
tions of the D vs. D(Mg) trends for olivine/liquid dis- 
covered by [2], The thermodynamic basis given by 
Jones [2] is seriously flawed because my analysis be- 
gan with the assumption that 

D = A/T + B (1) 

where D is the partition coefficient for some element 
between olivine and liquid, T is temperature, and A 
and B are constants. This faulty assumption can only 
be explained by Jones having had a very bad day. The 
correct starting point for any thermodynamic analysis 
should have been 

InD = A/T + B (2) 

Progress has been made in this area, however. 

The following discussions will distinguish Ni from 
other elements partitioning into olivine. 

Fe and Mn Partitioning. Linear parameteriza- 
tions of D(F e) and D(Mn) vs. D(Mg) can be simply 
explained because these regressions intersect the y-axis 
very near the origin. Therefore, the value of D(Fe) or 
D(Mn) can be approximated as zero when D(Mg) is 
zero. This is not the case for D(Ni), whose y-intercept 
is approximately -3.6. 

Differentiating Eq. (2) by 1/T gives 

01nD/0(l/T) = A (3) 

And in the Jones ([2] and Jones and Burnett [4] discus- 
sions of trace element partitioning, A is essentially 
equated with the heat of fusion of the mineral compo- 
nent of interest, i.e., fayalite, bunsenite, etc. 

Despite my several failed attempts to describe the 
D vs. D(Mg) trends, the thermodynamic analysis of 
those trends for Fe and Mn is quite simple. Consider a 
superposition of dlnD(Fe)/d(l/T) [for example] and 
dlnD(Mg)/d(l/T) giving 


01nD(Fe)/dlnD(Mg) = A(Fe)/A(Mg) (4) 

This can be expanded as 

[0D(Fe)/D(Fe)]/ [SD(Mg)/D(Mg)] = A(Fe)/A(Mg) (5) 
or 

dD(Fe)/dD(Mg) = K D (Fe/Mg) * A(Fe)/A(Mg) (6) 

A similar expression can be derived for D(Mn) vs. 
D(Mg). So the slopes of the Jones [2] regressions can 
today be calculated from known quantities: K D ’s for 
Fe/Mg and Mn/Mg and from heats of fusion (i.e., A’s). 
This calculation is only possible because the regres- 
sions for D(Fe) and D(Mn) vs. D(Mg) pass near the 
origin so that K D is nearly constant. 

Using heats of fusion for fayalite, bunsenite, and 
forsterite of 26.3 kcal/m, 32.5 kcal/m, and 27.2 
kcal/mole [3], resepctively, and K D (Fe/Mg) and 
K D (Mn/Mg) of 0.30 and .25, slopes on D(Fe) and 
D(Mn) vs. D(Mg) are calculated to be .29 and 0.30. 

The regression of experimental data for D(Fe) and 
D(Mn) give values of 0.30 and 0.26, respectively, in 
good agreement with prediction. Essentially, because 
the heats of fusion of fayalite, bunsenite, and forsterite 
are similar, the slopes of the D vs. D(Mg) regressions 
should nearly equal the K D of the exchange, as is ob- 
served. 

Ni Partitioning. The situation for elements such 
as Ni, whose D vs. D(Mg) regressions do not pass near 
the origin is more complicated. A consequence of eq. 
(6) is that the slope of a D(Ni) vs. D(Mg) diagram 
should not be constant because K D (Ni/Mg) is not con- 
stant. However, given that there are regions of D(Ni) 
vs. D(Mg) space that can be approximated as linear, I 
will use the approximation 

D(Ni) = a D(Mg) + b (7) 

to find an approximate expression for K D (Ni/Mg) 

K D (Ni/Mg) = a + b/D(Mg) (8) 

where a and b are published regression coefficients [3], 
so that eq. (6) becomes, after integration from 
D(Mg)=l to D(Mg): 

D(Ni) ={a (D(Mg) -1) + b lnD(Mg)} A(Ni)/A(Mg) (9) 


[91nD(F e)/9( 1 /T)]/[51nD(Mg)/S( 1/T)] = 



Figure 1 compares the published regression of [3] 
to the functionality of eq. (9). Except at large values of 
D(Ni), the comparison is fairly good, especially con- 
sidering that since eq. (8) is an approximation, eq. (9) 
must be as well. And, of course, no account has been 
taken for possible errors in the heats of fusion of the 
olivine components. Unlike the cases of D(Fe) and 
D(Mn), the A(i)/A(Mg) term is not approximately 
unity, but is closer to a factor of two [3]. 


Figure 1 



Can the functionality of eq. (9) fit real data? In 
Figure 2, 1 have regressed the combined data sets (L 3 ) 
of Longhi et al.[l] and Leeman and Lindstrom [5], 
using D(Mg) and lnD(Mg) as the independent vari- 
ables. The fit is quite good and is capable of explain- 
ing data at large values of D(Ni) that lie above the 
standard regression. 

Discussion. I have used simple thermodynamic 
arguments and data to explain linear relationships be- 
tween D(Fe), D(Mn), and D(Mg). The key to these 
relationships is the similarity of the heats of fusion of 
these olivine components. I have then shown that this 
same line of reasoning can be extended to D(Ni) vs. 
D(Mg) where differences in heats of fusion are large. 
However, this extension comes at the expense of a 
truly linear relationship. Nonetheless, an approximate 
functionality between D(Ni) and D(Mg) can be derived 
that serves to fit the experimental data. 

An immediate question is how/why a non-linear re- 
lationship between D(Ni) and D(Mg) has been viewed 
as linear for so long. A large part of the reason is that 
the non-linearity really only expresses itself at high 
values of D(Ni) where errors are often large. Also in 
the original Jones (1984) paper, there were data points 
that fell both above and below the regression line. 


With hindsight, though, it appears that the great major- 
ity of the points that fell below the regression came 
from a single study [6]. The reason for this is unclear. 
However, recent studies of D(Ni) for martian baalts [7, 
8] also fall on the multi-regression given here. How- 
ever, because these recent studies yielded intermediate 
values of D(Ni), they could equally well have been fit 
by the linear regression. Even so, it appears that the 
new regression is capable of describing D(Ni) in ter- 
restrial, lunar and martian basalts. 

It is also still unclear how D(Ni) will behave as 
D(Mg) approaches the orgin. However, taking eq. (9) 
at face value, its second derivative is positive at D(Mg) 
= 1 , implying that the function is concave down 
(“holds water”). This would also imply a flattening of 
the curve, relative to the standard linear regression. 
This flattening does not bode well for theories that 
postulate that Ni becomes incompatible in olivine at 
high temperatures. 

Finally, the heat of fusion of Co-olivine predicts 
that it will behave more like Fe and Mn than Ni, and 
this appears to be the case. The y-intercept of the 
D(Co) vs. D(Mg) regression is much closer to the ori- 
gin than that for D(Ni), although not quite as close as 
those for D(Fe) and D(Mn). 


Figure 2 
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